function G_cara = G_cara(V,gamma)
G_cara = -1/gamma*exp(-gamma*V);
end
